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ABSTRACT 

Recent observations have revealed that high surface-brightness, metal-rich debris is present over large regions 
of the Andromeda (M31) stellar halo. We use a set of numerical models to determine ■whether extended metal- 
rich debris is expected to exist in galaxy halos formed in a hierarchical ACDM universe. We identify tidal debris 
in the simulations according to the current best surface brightness detection limits in M3 1 and demonstrate that 
bright features in the outer halo regions usually originate from single satellites, with very little contamination 
from other sources due to the low probability of tidal streams from two overlapping accretion events. In 
contrast, high-surface brightness features in the inner halo often originate from multiple progenitors. We also 
compare the age and metallicity distribution of the debris with the well-mixed stellar halos in which they 
reside. We find that high-surface brightness tidal debris is produced almost exclusively by relatively high mass 
progenitors (M* 10^~^Mq) and thus is expected to be of moderate- to high-metallicity. Similarly, in our 
models the smooth inner halo is expected to be metal-rich as this region has been built up mainly from massive 
satellites. Our results imply that the stellar populations of substructure observed around external galaxies 
with current techniques should not resemble old and metal-poor dwarf spheroidal satellites, nor the underlying 
component of the stellar halo. 

Subject headings: galaxies: abundances — galaxies: evolution — cosmology: theory 



L INTRODUCTION 

One of the main goals of stellar population studies is to un- 
derstand the origin of the baryonic components of galaxies. 
Stellar populations contain a wealth of information - kine- 
matic, age and chemical abundance - that enables the recon- 
struction of star formation histories of the past and present- 
day galaxies. More recently, the growing interest in the 
topic of "near- field cosmology" (iFreeman & Bland-Hawthom 
I2OO2I) has raised the question of to what extent stellar popu- 
lations can be used to constrain cosmological models. In this 
paper we review recent stellar population data in M3 1 that ap- 
pear to challenge hierarchical models of structure formation 
and compare these data with results of numerical simulations. 

Galaxy halos are among the most interesting testing 
grounds for structure formation theories as they contain im- 
prints from various star formation episodes, starting with 
those that took place in the earliest formed structures. Pat- 
terns in the stellar populations in different regions of a galaxy 
halo, or between different halos, or even between halos and 
their neighboring satellite galaxies, can be analyzed and com- 
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pared with predictions of cosmological models. 

The old and metal-poor halo of our Galaxy has long 
been the p rototype f or models of g al axy formation 
(ISea rle & ZinnI 1 19781 iLaird et all 11988; Rv an & Norrii 
1I99II) . Thanks to advances in instrumentational tech- 
niques over the last couple of decades it has now become 
possible to resolve stellar populations in halos of other 
galaxies within our Local Group: Andromeda (M31) - 
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ellanic Clouds, M 33 and mo st of the dwarf sp he roidals 
dHarris & ZaritskyI 120011; iMcConacchie et al.l 120061 
Mighell&Rich ' 1996; Hurlev-Keller, Mateo & Nem3 
1998; Galla rt et al.i,1999i;,Harbeck et al.„200k.Dolphin.2002r 
Martin et al.l l2006h . Whereas the M31 and Milky Way halos 
exhibit similar stellar compositions (an underlying smooth, 
metal-poor component and a few metal-rich tidal structures), 
some local dwarf galaxies are quite different from the large 
halos. Some of these dwarfs have stellar populations that are 
on average younger and more metal-rich than the halo of our 
Galaxy, but overall display a great diversity in their stellar 
popul ations (Mateo 1998; Tolstoy et al. 2004; Dolphi n et al] 
l2005h . 

Recently, it has also become possible to analyze the stel- 
lar populations in some of the faintest galactic substructure, 
namely in the tidal debris. In contrast with the variety of pop- 
ulations displayed by neighboring dwarf galaxies, the streams 
discovered so far are striking in their uniformity. Unlike 
the predominantly old and metal-poor populations of nearby 
dwarf spheroidal satellite galaxies, the streams discovered 
so far are predominantly metal-rich with significant inter- 
mediate age components (Ferguson et al. 2005; Brown et al] 
l2006bi) . Moreover, the debris has distinct ages and metal- 
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licities even in comparison with the surrounding halo in 
which it resides. In the Milky Way, the Sagitta rius stream 
- the large st substru cture discovered in the halo jlbata et alj 
[1997; Maie wski et al . 1999, 2003) - has significantly younger 
stars than the stellar halo and an average metallicity 
that is a few dex higher ( Yannv et al.l l2000l: llvezic et al.l 
2000; Martmez-Delgado et al 
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Monaco et al...2006; .Belo kurov et al. '2006^ An even larger 
stream dominates the halo of M31 (Ibataet al. 2001). This 
so-called "giant southern stream" (GSS) has a mean metallic- 
ity of [Fe/ H1 = -0.5, higher by about 0.2 dex than the inne r 
M31 halo dFerguson et all l2002t lOuhathakurta et"an l2006l) . 
and by about 1 d ex than t he metal-poor M 3 1 halo component 
dChapman et al.l 2006; K alirai et al.ll2006l) . The stream also 
has a larger fraction of intermediate-age stars dBrown et al.l 
i006a). Other pieces of debris discovered in M31 in- 
clude the "NGC205 loop", the "Gl clump", the "North- 
East shelf", the "No rthern spur" dFerguson et al.ll2002ll200l 
iMcConacchie et all 12 004). etc. - the latest inventory done 
by 'Ib ata et all (l2007l) in M3 1 lists about a dozen or so tidal 
features. Detailed studies of different features in the color- 
magnitude diagrams in these regions suggest that the metal- 
licities are uniform and high, [Fe/H] ~ -0.4 dFerguson et al.l 
l200llFariaetani2007h . 

Are all the data consistent with the hierarchical model for 
structure formation? At first sight, the diversity of stellar pop- 
ulations in the halos, tidal streams and surviving satellites is 
suggestive of a stochastic assembly. However, some of the 
data have been given another interpretation. For example, it 
has been argued that the differences in chemical abundances 
between neighboring satellites and the Milky Way halo ex- 
clude thes e satellites from the list of possible halo progeni- 
tor types dVenn et al.ll2'004 ). At the same time, one can use 
the uniformity of chemical abundances in M3 1 's spheroid as 
an argument against the hierarchical scenario: if the halo 
formed from differen t pro genitors, wh^is the debris so uni- 
form in metallicity ( Ferguson et al.ll2005]) ? Other recent ob- 
servations claiming to challenge the hierarchical scenario use 
as arguments the apparent lack of very metal-poor stars in 
some of the surviving dwarf galaxies but not in the stellar 
halo (Helmi et al. 2006), and the metal-poor component of the 
M3 1 extending up to lar ge distanc es in an apparently homo- 
geneous fashion (Ibata et al.ll2007h . 

From the above arguments, it is clear that better theoret- 
ical expectations are needed in terms of the distribution of 
stellar populations in dwarf galaxies, stellar halos and in the 
debris. In a recent series of papers we have begun to com- 
pare the results of theoretical models with available data, us- 
ing as a basis a set of simulations of stellar halos formed in 
a ACDM cosmology. We found that the observed differences 
in [a/Fe] abundances between the surviving s atellites and the 
halo arise na turally in the ACDM framework (Rob ertson et al.l 
120051; iFonFet al. 2006a) . We have also shown how the phase- 
space structure of current debris can be u sed to recover the 
merger histories of AC DM stellar halos dFont et al.' '2006bl: 
iJohnston et al.llin~prep.l) or identify h alos whose me rger histo- 
ries are atypical for this cosmology dJohnston et aLllin prep.1) . 

In this paper we make a direct comparison between be- 
tween the stellar populations in our models and those in 
the halo of M31, down to a level of surface brightness 
similar to current detection limits. We choose M31 be- 
cause this galaxy has benefited in recent years from a sys- 
tematic study of its stellar populations, using either deep 
HST photometry, deep field spectroscopy or dedicated wide- 
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Gilbert et al. 2006, 2 007MKaUrai et al.l2006l;llbata et al 
Tanaka et al. 2007). As a result, some aspects of the 
stellar populations in the M3 1 halo are known better now than 
those in the Milky Way, where observations target mainly the 
Solar neighborhood. In addition, M3 1 exhibits a wide spread 
distribution of metal-rich debris throughout its halo. Our aim 
is to investigate whether this distribution is consistent with 
the hierarchical merger model. We use the age - metallicity 
relation (AMR) as an indicator for the formation of stars in 
our systems (both in the smooth stellar halo and in the tidal 
streams). We then compare the AMRs from our simulations 
with those obtained from the available stellar population data 
in M31. 

We also seek to constrain the origin of the bright, metal- 
rich debris in stellar halos (and in particular in M31). In ob- 
servations it is often difficult to disentangle multiple stream 
components in the debris. In some cases the stars may be- 
long to overlapping orbits of a single satellite (as it is now 
becoming evident in the case of the Sagittarius stream in the 
Milky Way). In others, the "debris" could be just the per- 
turbed or warped outer disk - as it may be the case for some of 
the re cent detections near the edges of M3 1 and Milky Way's 
disks (Momanv et all2006l:lFaria et al.l2007l) . When available, 
kinematic dat a can be u sed to separate the streams in phase- 
space ( Hel mi & de Ze euw 2000), and the modelling of orbits 
can further connect stellar features that are physically sepa- 
rate on the sky. Additionally, it is conceivable, although not 
yet proven, that stellar populations can be used as "tags" for 
the individual satellites from which they originate. In this pa- 
per, we use a statistical sample of bright tidal debris fields 
to provide a quantitative prediction for their origin. We esti- 
mate the average number of satellites per field with significant 
contributions to the bright debris, both in the inner and outer 
simulated halos. We also analyze the average metallicity of 
the bright debris and found that it is in good agreement with 
the observations in M3 1 . 

The paper is structured as follows: In section|2]we summa- 
rize the models . In section[3]we present th e AM Rs of the main 
galaxy halos 03.1b and of tidal streams ( ^3.2l l. In section [33] 
we investigate the origin of the highest surface brightness fea- 
tures ( ^3.3. lb . their typical ages and metallicities ( ^3. 3. 2b . In 
Section|4]we make detailed comparisons with ava ilable M 3 1 
data, for the halo (EB and the GSS (g42]and S2^- In 
Section|5]we conclude. 

2. MODELS 

The set of stellar galaxy halos used here has been con- 
structed by iBullock & JohnstonI d2005l) and it comprises 1 1 
dark matter halos of mass Mi,,> = 1.4 x lO'^M© formed in a 
ACDM Universe (n„ = 0.3, 17a = 0.7, h = 0.7, and erg = 0.9). 
For our discussion we will also refer to these as L» galaxy 
halos. 

Full details about the set-up of the simulations can be 
found in the previous pa pers (Bullock & Johnston 2005} 
[Robertson e"taLll2005l;lFont et al.ll200 6a). Here we just briefly 
summarize the main characteristics of the models. The merger 
history of t he galaxies is based on an Extended Press-Schecter 
formalism dLacev & Cold [T993h and the evolution of each 
satellite is followed with a numerical simulation after accre- 
tion onto the main halo. The baryonic components are mod- 
eled with semi-analytical prescriptions. The infall of cold 
gas into each satellite is followed from the epoch of reion- 
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ization to the time of its accretion onto the main halo. The 
gas mass is translated into a stellar mass through a simple 
star-formation prescription and the stellar profile is assumed 
to follow a King (1962) model. Star formation ceases soon 
after accretion onto the main halo when the gas is assumed to 
be lost due to ram-pressure stripping. The chemical evolution 
of the accreted substructure is followed with the method of 
[Robertson et alj (i2005i) . which takes into account the enrich- 
ment from both Type II and Type la supernovae and models 
feedback as metal-enriched winds from both supernovae and 
intermediate mass stars. The models evolve with redshift, but 
are calibrated only at redshift z = to match the stellar mass 
- metallicity relation, -Z, and the stellar mass - circular 
velocity relation, M^-Vcirc, for Local Group dwarfs. 

When evolved to the present day, the models match the 
main physical properties of the two main galaxies in the Lo- 
cal Group and of their progenitors: the number of bright (> 
IO^Lq) surviving satellites per main galaxy, and the masses, 
density profiles and total luminosities of the stellar halos are 
similar to those found in the Milky Way and M3 1 . The aver- 
age metallicities of our simulated L*-type galaxy halos range 
between -L3 and -0.9. The model also retrieves the differ- 
ences in the [a/Fe] abundances between stellar halos and the 
surviving satellites. More details about the spatial distribution 
of chemical abundances ([Fe/H] and [a/Fe]) in these stellar 
halos can be found in Font et al. ( 2006a 

3. RESULTS 

3.1. Predicted properties of stellar halos 

The AMR of a stellar population is a useful tool for gaining 
insight into the chemical evolution of the system. In a closed 
box system, it is expected that the oldest populations will 
be the most metal poor whereas the youngest ones the most 
metal rich. The slope of the AMR can vary, being generally 
steeper for the initial stages of the evolution. The AMR will 
be monot onic only for is olated sy stems that form ed mono- 
lithically (iPagel & Patche tt 1975; B inney & Merri field 1998|, 
while systems that formed by accretion of satellites will dis- 
play a composite functional form, with multiple metallicity 
components for the same age. In this section we use our mod- 
els to investigate to what extent the AMR of a stellar halo can 
be used to place constraints on the galaxy mass assembly his- 
tory. 

Figure [T] shows the AMR of our "typical" stellar halo, 
constructed by averaging over all 11 halos in our sample. 
Each symbol represents a cell in the 2-dimensional age - 
[Fe/H] grid where there are any stars' \ with the area of the 
symbol corresponding to the mass fraction contributed by the 
stellar population to the total stellar mass of the averaged halo. 
Only the inner 50 kpc of the stellar halos are considered for 
this plot. 

The AMR in Figure [T] has a number of distinctive charac- 
teristics: 1) the AMR has a single dominant sequence rather 
than being composed of multiple ones; 2) the sequence is 
representative of satellites that are massive, i.e that enrich to 
[Fe/H] > -1 values in just a few Gyrs; and 3) the ages of the 
stars are typically old, > 8 Gyr. These characteristics are 
a direct consequence of the merger histories of inner halos 
assumed for our ACDM models, and w hich have been de- 
scribed in more d etail in previous papers dBuUock & JohnstonI 
l2005t iFont et alJi2006a.b) . The dominant AMR sequence oc- 

' ' We henceforth consider all stars in a given grid age - [Fe/H] cell to be a 
distinct stellar population. 
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Fig . 1 . — AMR for the average of all 1 1 halos in our sample. Only the inner 
50 kpc regions of the halos have been considered here. The area of symbols 
varies with the mass fraction of stellar populations to the total mass in the 
field. (No population contributes, by mass, more than 25% to the halo.) 
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Fig. 2. — AMR for the outer region (50 < R kpc) of the averaged simulated 
halo. Symbol areas coiTespond to the stellar mass fraction contributed by the 
population, as in Figure [l] The fraction of each population is normalized to 
the stellar mass of the entire outer halo region. 



curs because the bulk of the inner halo is formed from a 
few massive satellites (10^ - IO^^Mq), even though the main 
halo accreted satellites with a broader range of stellar masses 
(M, - IO'^-IO^^Mq). In addition, massive satellites sink more 
rapidly into the center due to a more effective drag of dynam- 
ical friction. The major contribution to the total stellar mass 
of massive satellites also explains the clear signature of rapid 
metal enrichment (more details about the AMRs of individual 
satellites are given below). Few stars have ages < 8 Gyrs, and 
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Fig. 3. — AMRs of individual satellites of various masses and accreted at different times. The low mass satellites chosen here have M, ~ 10* Mq, while the 
high mass ones have M* of a few IO^Mq. The accretion times are of 1—2 Gyr for the early events and 6-7 Gyr for the late ones. Symbol areas correspond to 
the stellar fraction contributed by the population, as in Figure [T] The stellar mass fractions are normalized to the total stellar mass of each satellite. 



this is because early events dominate the inner halo accretion 
history. 

Spatial variations in the stellar populations properties 
across hierarchically-formed halo are also of observational in- 
terest. Figure |2] shows the AMR of the outer regions of halos 
(R > 50 kpc), using the stellar population data of the same 
averaged halo as in Figure [1] 

The larger spread in ages in the outer halo is a consequence 
of th e late assembly of this region compared wi th the inner 
one ( Bullock & Johnstonll2005FlFont et al.ll2006al) . Given the 
evolution with time of the typical mass of the accreted sys- 
tems, the outer halo is formed from satellites with higher 
masses, ~ 10^- IQ'^Mq, than the inner halo. Of note, 
is the presence of very massive satellites in the outer halo, 
M, ~ IO^Mq, which extend the outer halo AMR towards 
[Fe/H] ~ values. 

Our results indicate that an averaged sample of the 
outer halo should display a significantly larger fraction of 
intermediate-age versus old stars than the inner halo. The 
metallicity distribution function should also have a broader 
peak in the outer halo. 

3.2. Predicted properties of individual satellites and streams 

By construction our simulated satellites have a monotonic 
chemical enrichment. This limitation results from a num- 
ber of simplifying assumptions - for example, the merger 
history of a satellite is approximated with a general formula 



(IWechsler et al.l2002h and not followed stochastically, and the 
gas and metals are assumed to be instantaneously mixed af- 
ter each star formation episode. The predictive power of our 
models for individual satellites (and streams) is therefore re- 
stricted. The individual AMRs will underestimate the scatter 
in [Fe/H] for stellar populations of the same age. However our 
models are fairly reliable in matching the mean metallicity of 
satellites and the overall AMR sequence. 

Figure |3] shows four cases of individual satellites in the 
range of stellar masses and accretion times that are typically 
found in our simulations. Despite the limitations in modelling 
the entire scatter, some trends are obvious. For example, satel- 
lites have a smaller [Fe/H] spread at a fixed age than the outer 
halo but the spread is comparable to that of the inner halo. 
Also, the overall shape and maximum [Fe/H] of the inner halo 
AMR are similar to those of massive satellites (M* > IQ^Mq). 

The AMR of the outer halo contains individual AMR se- 
quences that are more clearly separated (Figure |2|i than those 
in the inner region of the halo (Figure [T]i. As each AMR se- 
quence corresponds to a single satellite or a group of satellites 
with similar properties, this suggests that using AMRs as a 
tool to disentangle the accretion histories of stellar halos may 
be more feasible in the outer rather than in the inner regions 
of the halos. 

3.3. The Extended, High Surface Brightness Debris 
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Fig. 5. — The origin of liigli surface brightness fields (/iv < 30 
mag/arcsec^): a single source (case A), two dominant sources (case B), mul- 
tiple less dominant sources (case C) and diffuse material originating in many 
minor sources (case D). Top panel: All bright fields across the halo are se- 
lected; Bottom panel: Only the bright fields in the inner R < 10 kpc of the 
halo are selected. 



Fig. 4. — Top: A stellar halo mapped in surface brightness down to 34 
mag/arcsec^ (this is halo H2 in Font et al. (2006b)). The region shown here 
extends 150 kpc on the side. Middle and bottom: Four examples of kinemat- 
ically cold streams surviving at present time in the halo above. The sfi'eams 
have been shrunk to a quarter of their size, for illustration purposes. 



3.3.1. Origin 

What is the origin of the bright, metal-rich debris extend- 
ing over large areas in the M31 halo? Was a single satellite 
responsible for polluting the M3 1 halo with metal-rich stars 
or were they several? Here we use our models to constrain 
the formation history of this debris, assuming that it formed 
entirely from merging events. 

The upper panel in Figure |4] shows the surface brightness 
map of a stellar halo with M31-like properties. The map 
shows only the surface brightness down to the limit of /iy = 34 
mag/arcsec^, which is roughly the deep est limit achieved cur- 
rently by wide fie l d surveys of the M 3 1 (iFerguson et al.ll2002t 
llrwinet all [20051; llbataet al.ll2007l) . The bottom panels in 
Figure |4] show the surface brightness maps of the main four 
kinematically cold streams surviving at the present time in 
the simulated halo. 

Down to the chosen surface brightness cut, the spatial dis- 
tribution of the simulated debris is similar to what has been 
found so far in M31 and in the Milky Way: there are a 
few bright kinematically cold streams (the counterparts of the 
Sagittarius stream, the GSS or the NGC205 loop), the rest of 
the debris being in the form of much fainter streams or diffuse 
overdensities. 



We define as "high surface brightness" all debris that is as 
bright (or brighter) than the current detections of su bstruc- 
ture - iJLv ^ 30 mag/arcsec^ for the GSS in M31 (Ibat a et al] 
1200 lb . We therefore select for our analysis all fields which 
have < 30 mag/arcsec^ and that are located on the pro- 
jected area of a halo along the line of sight. These mock fields 
have the same area as the pixel size area of the halo map in 
Figure m i.e. 0.4 x 0.4 kpc. The depth of the field along the 
line of sight is of the same order as the virial radius of the 
main galaxy, i.e. ~ 300 kpc. The total surface brightness in 
each pixel is obtained by integrating over all populations in 
the field. 

We then divide all bright fields into four disjoint cases, 
based on the formation history of the stellar overdensity in 
that field: 

- case A: The overdensity originates mainly in a single 
satellite. The condition we impose is that the satellite con- 
tributes a mass fraction / > 80% of the total stellar mass in 
the field; 

- case B: The overdensity is formed from contributions of 
two main satellites. The condition is that these satellites com- 
bined contribute more than 80% of the stellar mass in the field; 

- case C: The overdensity is formed from one or more satel- 
lites with less significant contribution(s). Specifically, we re- 
quire that there is at least one satellite contributing 30% < f < 
40%, but no satellite contributes more than 40% to the total 
stellar mass in the field; 

- case D: The overdensity is formed from a multitude of 
minor satellite contributions. No satellite is contributing more 
than 30% of the total mass in the field. 

The top panel of Figure|5]shows the frequency of each case 
among the bright fields over the entire area of the stellar halo. 
Cases in which the overdensity originates in a single satellite 
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(A) occur about 70% of the time, while cases B, C and D oc- 
cur less frequently, each at about 10% of the time. The actual 
percentages vary slightly from halo to halo and depend on the 
imposed surface brightness limit - for example, the frequency 
of case A increases if we select an even brighter debris. How- 
ever, case A consistently dominates across the halos. The dis- 
tribution of these cases begins to reverse in the inner regions 
of the halos, where phase-mixing is more important. The bot- 
tom panel in Figure |5] shows the same information as the top 
panel, but for the inner 10 kpc of the halo. Here, the dominat- 
ing cases are those in which the debris originates in three or 
more satellites. 

These results have a number of implications for the inter- 
pretation of high surface brightness overdensities in galaxies. 
For example, if the overdensity is relatively distant from the 
center of the halo, R> 10-20 kpc, it is likely to originate in a 
single satellite, with little contamination from other satellites 
in the field. This result holds regardless of whether the feature 
is kinematically cold or if it is diffused in phase-space. For the 
inner halo, only a third of the bright features are created by a 
single satellite, the majority of them originating from super- 
position of streams from three or more satellites (cases C and 
D). 

3.3.2. MetalUcity 

The existence of the metallicity - luminosity rela- 
tion (Garnett & Shields 1987; Skillman, Kennicutt, & Hodee 
1989;.Brodie & Huchra..l991;,Zaritskv. Kennicutt. & Huchra 
19941) and the corresponding metallicity - mass relation 
(iLequeux et aP 1 1 9791: llYemonti et all l2004t lErb et all l2006l) 
for galaxies indicates a connection between the nature of the 
satellite progenitors and their metallicities. We might expect a 
corresponding connection to exist between the surface bright- 
ness of the current debris and its overall metallicity. Given 
this expectation, it is useful to examine whether high surface 
brightness debris in our models (as analogue of the M3 1 ob- 
servations) indeed originates in satellites of a typical stellar 
mass, and if so, what is the predicted average metallicity of 
the debris. 

In our models, satellites relatively massive, 
and accreted 10-12 Gyr ago result in sufficientl y bright de- 
bris to b e detected, ^.y ^ 30-32 mag/arcsec^ (iJohnston et al.l 
lin prep.D . As the surface brightness of the debris depends on 
both the mass of the progenitor, and on the time of accre- 
tion, satellites with ^ lO^'^M© accreted less than 10 Gyr 
ago are much brighter, ^ 2%- 30 mag/arcsec^. Satellites 
with stellar masses greater than lO^M© assembly later and 
are therefore accreted later (< 8 Gyr ago), resulting in debris 
with the highest surface brightness, /iy ~ 26-28 mag/arcsec^ 
(iJohnston et al.l lin prep.V Therefore, the average metallici- 
ties of high surface brightness fields originating from a sin- 
gle satellite are expected to range between -1 .4 and (which 
are average values for satellites with stellar masses of 10^ and 
1O^M0, respectively; values are quoted at the time of accre- 
tion). 

A high surface brightness in a given field can be achieved 
not only from the contribution of a single relatively massive 
satellite, but also from the superposition of debris originating 
in several massive satellites (eg. if 80% of field stars originat- 
ing in a 1O^~'^M0 satellite reach the detection limit, i.e a case 
A, then a combination of two 40% fractions originating in two 
separate 1O^~**M0 satellites are equally bright, i.e a case B). 
In the case of superposition of massive satellites, the average 
metallicity of the high surface brightness field is similar to the 



average metallicity of the individual satellites, i.e. the fields 
are expected to be metal-rich. 

As we will discuss below, most of the observed debris in 
M31 is metal-rich. The average [Fe/HI i n the GSS is - -0.5 
(Ferguson et al. 2002'; 'Br own et al.l '2006b) and other debris 
throughout the M31 halo have similarly hi gh metallicities, 
rFe/HI ^ - -0.4 (iFerguson et allllOOl iF^iaet al. 2007). Ad- 
ditionally, |GifcerteialI^2007b show that the M31 "spheroid" 
fields and several fields along the minor axis are consistently 
metal-rich, with and without the cold GSS component. The 
lesson learnt from our models is that the metal-rich debris cur- 
rently observed in the inner halo of M3 1 may not have a sin- 
gle common progenitor, and moreover, each individual inner 
halo field may contain metal-rich components from multiple 
sources. 

4. APPLICATIONS TO OBSERVATIONS 
4.1. TheM31 Halo 

iBrown et al.l (l2006al[bl l2007h have derived the star forma- 
tion histories at various locations in the outskirts of M3 1 . The 
resulting ages and metallicities in these fields can be com- 
pared directly to our model predictions. A concern is that 
some of Brown's fields lie along the minor axis of M31, a 
region that has been shown to cont ain stars from m ultiple 
sources, including the disk and GSS (llbata et al.ll2007h . The 
most recent field analyzed lies at 21 kpc along the minor axis 
and is expected to have an insignificant contribution from thin 
disk stars o r from the giant extended disk structure detected by 
llbata et all (l2005i) (for e vidence again s t cont amination from 
the latter component, see iGilbert et al.l (l2007l) ; also, note that 
the extended disk component may have an accretion origin 
as well (Penarrubia, McConnachie & Babul 2006)). On the 
other hand, the 2 1 kpc halo field may still be contaminated by 
sti-eam debris (iGilbert et a l.. 2007 ). 

We note, however, that from the point of view of hierarchi- 
cal models one cannot easily differentiate between a "true" 
or "clean" halo and a halo contaminated by debris. In this 
context, stellar halos are expected to form through the super- 
position of tidal debris and show various degrees of phase- 
mixing. Therefore the only concern in our comparison is not 
to mistake a field that is dominated by a given stream for one 
that is representative of the whole halo. The cautionary re- 
m ark regarding our further comparison with the 21 kpc field 
of iBrown et al.l ((2007) is that the presence of a cold stream in 
the inner halo field appears indeed to be atypical in our sim- 
ulations. Whether this observed field is representative of the 
whole inner M31 halo, or of a typical Milky Way- type inner 
halo, remains to be decided by future observations. 

Our strategy for further comparisons with the observa- 
tions is to present averaged results rather than selected, well- 
matched fields. In this respect, our results do not emphasize 
the dispersion in stellar population data that can exist among 
individual halo fields. However, these differences are not ex- 
pected to be as l arge in the inner h alo as in the less phase- 
mixed, outer one (iFont et al.ll2006bl) . 

4.1.1. The 2 1 kpc Halo Field 

Figure m shows a com parison between the 21 kpc halo field 
of iBrown et al.l (12007 ) and our averaged inner halo region 
(bottom and top left panels, respectively). In the right pan- 
els we also show the AMR of two individual halos from our 
sample that have mass accretion histories that bracket the av- 
eraged one (halo H7 has a protracted mass accretion history. 
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whereas halo H6 has a short one. Halo H7 is also 0.5 dex 
more metal-rich than halo H6). 

The dominant presence of old stars in both observations 
and simulations is not surprising, as inner halos are expected 
to form within the first fe w Gyr from the assembly history 
(iBullock & JohnstonI 12005 ). The existence of intermediate- 
age stars in the inner halos is a new finding that can pro- 
vide important cons traints on the models of halo formation. 
iBrown et all (l2006bl) estimate that ^ 40% of the stars in the 
field 1 1 kpc along the minor axis in M3 1 are less than 10 Gyr 
old, while the mean age in that field is 9.7 Gyr. Brown et al.l 
(|2()07) find that the field at 21 kpc along the minor axis has 
significantly less stars with ages less than 10 Gyr (the fraction 
is about half of that in the inner field), and virtually no stars 
with ages less than 8 Gyr The distribution of populations in 
the latter field seems to be more into agreement with our mod- 
els, which predict that about 20% of populations in the inner 
halo have ages less than 10 Gyr. 

The AMRs of the simulated halos show non-negligible pop- 
ulations with intermediate metallicities -2 <[Fe/H] < -1 and 
ages younger than 10-11 Gyr. Although these populations 
represent only a minor contribution (< 1%, by mass) to the 
entire stellar halo, they can provide important information 
about the halo formation history. In our models, stars with 
these ages and metallicities predominantly come from lower- 
mass, metal-poor halo progenitors, and hence are a probe of 



the accretion of dwarf galaxies at the lower end of the mass 
function. (Note: the lower metallicity of low mass halos fol- 
lows because we i mposed that dwarfs accreted today follo w 
the M^- Z relation ^ekel & SiUjT986l:lDekel & Wool2003h .) 

The old and metal-poor stars with ages > 10 Gyr and 
[Fe/H] < -2, located in the lower right corners of the AMRs, 
are good tracers of the early stages of star formation in galaxy 
halos or in their progenitors. Both our models and the ob- 
servations show non-negligib le fractions of such metal-poor 
stars (note: the i Brown et al.l (2007) AMR shows only stars 
with [Fe/H] > -2, but this may only be an artifact of the lim- 
ited choice of metallicities in the isochrone dataset). At the 
moment there is a great interest in finding the most metal- 
poor stars in halos of galaxies and in understanding their ori- 
gin. Apart from the evidence that the M31 halo has some 
very metal-poor stars dBrown et all 12004 ', 2006b, 2007), the 
search es also target the Milky Way halo (Beers & Christlie3 
l2005h . In our models most of the stellar populations form in 
small sub-halos at high redshift and enter the main halo in 
relatively massive satellites. Improved models that track the 
early episodes of star formation and metal enrichment in more 
detail will be able to provide better constraints on the origin 
of this type of populations (see Tumlinso n et al, in prep.). 

Interestingly, the observed AMR of iBrown et aP (l2007h 
contains some stars that are more metal-rich than the stars 
of the same age in our simulated halos (see, in particular, the 
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[8 Gyr,12 Gyr] age interval). This rapid enrichment is indica- 
tive of stars originating in a fairly massive satellite, of 10^ 
or a few 10^ Mq in stellar mass (see, for example the AMRs 
of massive satellites in Figure O. The reason that our sim- 
ulated halos do not show such signatures is that most of the 
progenitors of these well-mixed regions have stellar masses 
of IO^'^^Mq or less. According to our models, the accretion 
of a ^ 10^"^ Mq satellite is a rare event it it were to occur 
more than 8 Gyr ago. A more recent accretion of such a mas- 
sive satellite is possible, however the chances that the satellite 
is fully disrupted such that it mixes with the rest of the halo are 
small, unless its orbit is extremely radial. A possible explana- 
tion for the origin of these old and metal-rich populations is 
that they ori ginate in the GSS (recall that the 21 kpc field of 
iBrown et alj ( 1200711) is contaminated by the GSS; and that the 
progenitor of the GSS is believed to have been fairly massive 
and on a radial orbit. See also the discussion of the GSS data 
in sections ^4X11 and ^4Jj\ . 

There are also a number of discrepancies be tween our mod- 
els and the observations of lBrown et al.l (12007!) . The stars with 
ages < 5 Gyr and [Fe/H] > -1 in the top left corner of the ob- 
served AMR cannot be matched by our models. It is, however, 
unlikely that these extremely young and metal-rich stars be- 
long to the halo. We hypothesize that they are either caused 
by a very minor contamination by the disk in this field, or that 
they belong to a recently accreted stream. The reason that our 
simulated AMRs do not have this young and metal-rich com- 
ponent is because, on average, streams of satellites accreted 
less than 6 Gyr ago are not found in this inner regions of our 
simulated halos. 

Other mo r e extreme ty pes of populations inferred by 
iBrown et alj ( l2006bl 120071) in the halo of M31 are also not 
easily retrieved in our models - for example, the trace popula- 
tions of young and extremely metal-poor stars (ages < 10 Gyr 
and [Fe/H] < -2) in the bottom left corner of the observed 
AMR or the old and metal rich stars (ages > 10 Gyr and 
[Fe/H] > -0.5) in the top right corner of the same diagram. 
The reason for this mismatch may be that the modelling of 
the gas and stellar physics in our models is too simplistic to 
account for these populations. However it is difficult to under- 
stand what type of satellite systems would have such unusual 
paths of chemical evolution. In particular, if the young stars 
with metallicities [Fe/H] < -2 do indeed have a merger ori- 
gin, it implies that the progenitor satellites either formed only 
recently and in relative pristine environments, or that the mix- 
ing of metals inside these satellites was very inefficient. We 
note, however, that while the presence of these ex treme pop- 
ulations was derived from the iBrown et alJ (l2006b) analysis, 
the authors note themselves that they may not actually be re- 
quired. Indeed, IBrown et al.l (12006 b') showed that old metal 
rich stars could easily be excluded from their model without 
diminishing the goodness of fit. Furthermore, they point out 
that what appears as a young metal poor population in their 
CMD could at least partially be explained as a population of 
blue stragglers. 

4.1.2. The Outer Halo 

Figure |2] also shows that the metal-poor populations 
([Fe/H] < -2) extend out to large distances in the halo, in 
agreement with recent observations in M3 1 which detect an 
underlying metal-poor halo component out to ^ 150 kpc 
jChapman et al.!)2006t [Gilbert et al.ll2006l: iKahrai et al.ll2006l: 
llbataetal.2007fK 

We find that the ratio of metal-poor stars to metal-rich stars 



is higher in the outer halo than in the inner one'^. This re- 
sult appears to be in agreement with observations that find 
that the overall me tallicity of stars decrease s with radius in 
the MSl's spheroid dBellazzini et al .' 2003; Ka lirai et alj2006l; 
Brown et al. 2007). In contrast, the density profiles derived by 
Irwin et al.. (2005); Ibata et al. (2007) for the metal-poor and 
metal-rich populations in M3 1 appear very similar, suggest- 
ing that the ratio of these populations is fairly constant from 
inner to the outer halo. 

We also predict that the outer regions of the halos should 
have, on average, younger ages than the inner ones, as a con- 
sequence of the in side-out growth of these halos. The two 
fields examined by lBrown et alJ (l2007h at 11 kpc and 21 kpc 
along the minor axis of M3 1 show an inverse trend. However, 
it is likely that the excess of younger stars detected in their 
inner field is the result of contamination by the disk, rather 
than a genuine trend (however, note that Brown and collabo- 
rators argue against disk contamination in the 1 1 kpc field). It 
is certainly desirable to have a detailed analysis of both ages 
and metallicities in more fields across the M31 halo before 
drawing a definitive conclusion about the age and metallicity 
gradients in this system. 

4.2. The Giant Southern Stream 
4.2.1. Stellar Populations 

Figure |7] s hows t he AMR of the GSS field derived by 
iBrown et alj (l2006bh alongside the AMRs of several simu- 
lated streams. We have chosen for this comparison only 
streams that have massive progenitors, M„ > 10 Mq, as sev- 
eral observational and theoretic al studies indicate that this is 
the cas e of the G SS progenitor (Ibata et al. 2004; B rown et alj 
,2006ai;lFont et aU2006c) . For example, we find that only pro- 
genitors with total stellar masses greater than ~ 5 • IO^Mq 
are able to match the most metal-rich portion of the observed 
AMR. Among these massive satellites, we find that those ac- 
creted ~ 6 - 7 Gyr ago provide a better fit to the observed 
AMR (eg. the bottom left panel in Figure |7]) and also result 
in tidal debris of comp arable surface bri ghtness as the GSS, 
liv ^ 30 mag/arcsec^ ("I bata et ani2001h . (A more detailed 
discussio n abou t constraining the time of accretion will be 
given in ^4.2.2b . 

The AMRs of satellites accreted 6-7 Gyr ago also have 
about 70% of their stars younger thar i 10 Gyr, in agreern ent 
with what has been found in the GSS dBrown et alj|2006bl) . 

Not all stellar population data can be matched by our mod- 
els. The metal-poor populations ([Fe/H] < -2) are not as sig- 
nificant in the simu lated AMR than in the stream field of 
iBrown et alj (l2006bl) . Also, as in the case of main halos, the 
most extreme populations detected in the GSS field, such as 
the top right, to p left and bottom left corners of the age - 
[Fe/H] grid of Brow n et"aLl (l2006bl) are inconsistent with our 
chemical evolution models. 

As discussed before, the mismatch between the observa- 
tions and the models can be attributed, in part, to the limi- 
tations of our current modelling of individual streams. Ad- 
ditionally, the observational data in the GSS field may con- 
tain stars from other components, such as the dis k, other tidal 
debris, or the background halo. For instance, iBrown et alj 
( l2006bl) estimate that about 25% of the stars in the stream 

This result does not imply that the absolute number of metal-poor stars 
in outer halo should be higher than in the inner halo. The density of metal- 
poor stosjsexpected to decrease with distance from the center of the galaxy 
(see lFont et sd\ )2006bD ). 
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field could belong to the underlying halo. In addition, as pre- 
viously discussed in ^4.1.11 some of these extreme popula- 
tions, such as the old and metal-rich ones, are not necessarily 
required for a good fit to the observations (although the young 
and metal-poor ones seem to be required). 

4.2.2. Constraining the Time of Accretion 

Although the constraints on the mass of the GSS progenitor 
are fairly good, the time of accretion is largely uncertain'^. 
Here we discuss the possibility of putting some constraints on 
the latter by using our star formation history models. 

Even though stellar population data contain direct infor- 
mation only about the star formation history of the progen- 
itor galaxy, indirectly, they can also provide some constraints 
on the external mechanisms that can shut off the star for- 
mation. The accretion of the satellite onto the main halo is 
such a mechanism that may affect the star formation process 
within the sa tellite, slowing it down or stopping it completely. 
iBrown et alj (l2006b) suggested that the cessation of star for- 
mation in the GSS ~ 6-7 Gyr ago might have coincided with 
the accretion of the progenitor satellite onto the M3 1 . For ex- 
ample, if the satellite were accreted earlier than 7 Gyr ago, 
ram-pressure stripping along this highly eccentric orbit could 
have removed the gas reservoir and therefore stopped the star 

The time of accretion is defined here as the time when the satellite enters 
the virial radius of the host galaxy. 



formation earlier than it is observed. Similarly, if the satel- 
lite were accreted more recently, star formation would have 
continued several Gyrs after what is observed, unless some 
intrinsic mechanism would have stopped it before the accre- 
tion. 

However, a time of accretion of 6 - 7 Gyr ago is at odds 
with current numerical modelling of the satellite orbit that 
fits the kinematic data. The simulations find that the pro- 
genitor satellite has spent less than 1 Gyr orbiting from the 
farthest location along the stream until the current pericen- 
ter (Ibataetal. 2004; Fardal et al. 2006; Fontetal. 2006^ 
iGeehan et a l. 2006; F ardal et al.l2007l) . and assuming that far- 
thest edge of the stream is the point where the satellite was 
accreted, it implies that the satellite must have been accreted 
within the last ~ 1 Gyr. 

Although the stellar population data cannot provide strin- 
gent constraints on the dynamical evolution of the satellite, 
the arguments above point towards a possible conflict with 
the numerical estimates. We examine below two possibilities: 
first, that the progenitor was accreted recently (< 1 - 2 Gyr 
ago), but the star formation stopped several Gyr prior to that; 
and second, that the satellite was accreted 6-7 Gyr. 

In the first case, the question is how likely is that star forma- 
tion stops several Gyrs prior to the accretion, particularly in a 
fairly massive satellite as it is inferred to be the case for the 
progenitor of the GSS. In Figure [8] we investigate this sce- 
nario quantitatively, by selecting all satellites accreted less 
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than 7 Gyr ago and plotting the fraction of stars formed in 
the 4 Gyr time interval before accretion. The satellites are 
divided in three categories depending on their time of accre- 
tion: satellites accreted about 5 -7 Gyr ago, satellites accreted 
recently (< 2.5 Gyr), and satellites accreted at intermediate 
times. The fraction of stars is expressed function of the to- 
tal stellar mass of the satellite at the time of accretion. For 
more massive satellites, M» > IO^Mq, - as it is assumed for 
the progenitor of the GSS - at least 40% of all their stars are 
formed within the 4 Gyr before accretion. The implication 
is that, if the progenitor were indeed accreted within the last 
1-2 Gyr, the observations should detect a ^ 40% fraction of 
young stars (ages < 6 Gyr). Current observatio ns, however, 
do not detect such a high fraction of young stars (iBrown et al.l 
|2006b). 

Therefore, the results based on our models seem to concur 
that, for a massive satellite, a time of accretion of 1 - 2 Gyr ago 
is difficult to reconcile with a cessation of star formation 4 Gyr 
prior to that. The caveat, of course, is that the star formation 
prescription in our models may be too simplistic. The alter- 
native is that the progenitor of the GSS was indeed accreted 
6-7 Gyr ago. In this case, it is likely that the satellite orbited 
more than once in the M3 1 halo, as the time taken to complete 
an orbit with the current kinematic properties is much less 
than 6 -7 Gyr (see the numerical simulations of llbata et al 



(|2004V'Fardal et al.' (2006|): lFont et all (l2006cl) : lGeehan et al 



12006); Fardal et al. (2007|)). Also, it is also unlikely that a 
tidal stream that detached 7 Gyr ago would still survive to- 
day as a coherent feature. Therefore, if the satellite was ac- 
creted 6-7 Gyr ago, the stars in the current stream must have 
been tidally stripped from the parent halo relatively recently. 
It is possible that a resolution between the results of numer- 
ical modelling and stellar population data can be achieved if 
the starting point of simulations is moved beyond the edge 
of the currently detected stream. Extrapolating the numerical 



parameter space in regions which are not constrained by the 
data would be difficult, but certainly worth doing. The task of 
constraining the orbital parameters could be greatly eased if 
future observations find a continuation of the stream beyond 
the current edge, either along this orbit or on a previous one. 

Our current simulations can provide some predictions on 
how faint the streams of satellites that undergo multiple orbits 
should be. The small panels at the bottom of the Figure |4] il- 
lustrate the surface brightness of a few streams formed when 
the parent satellite is at its first and then second passage inside 
the main halo. The streams are brightest when they just form, 
then faint with time. Therefore the streams that formed in a re- 
cent passage of a satellite inside the main halo are expected to 
be brighter than the streams that survive from previous a orbit 
of that satellite (however, the effect is attenuated by the fact 
that more stellar material is stripped at earlier passages of the 
satellite, therefore the initial surface brightness of the streams 
is higher for the streams formed earlier). We find that while 
the brighter (i.e. more recently detached) streams have an av- 
erage surface brightness of 25-26 mag/arcsec^, the fainter 
counterparts have, on average, fiy ^ 26-27 mag/arcsec^. 
These surface brightness levels should be achievable with cur- 
rent detection limits. 

5. CONCLUSIONS 

In this paper we used a set of simulations of stellar halos 
of galaxies to study the composition of their stellar popu- 
lations and provide comparisons with the observations. The 
test case for our comparisons was M31, given the wealth of 
data that exists for this galaxy. However, we emphasize that 
the theoretical results derived here are applicable to any large 
^ L* galaxy and therefore provide a framework for interpret- 
ing the substructure of such galaxies that will be uncovered in 
upcoming surveys. 

The main conclusions of our study are as following: 

- The global stellar populations of halos assembled hierar- 
chically will vary according to their different accretion his- 
tories. For example, the presence of metal-rich debris in a 
galaxy like M3 1 suggests a history that included the accretion 
of at least one massive satellite. Conversely, the more metal 
poor halo of the Milky Way is evidence that its progenitor 
satellites were less massive than those of M3 1 . 

- The homogeneous high m etallicity in the inner M31 halo 
found by lBrown et al.l (l2007h can be naturally explained in the 
context of hierarchical merging. According to our models the 
homogeneity of stellar populations occurs because the inner 
halo is dominated by a few massive satellites, ~ 1O^~**M0 , 
which had very similar star formation and chemical evolution 
histories. The predicted imprint of these satellites in the inner 
stellar halo is the presence of old stars with [Fe/H] ~ -1 0. 

- The outer halos in our models tend to have a larger spread 
in the ages and metallicities of stellar populations than the 
inner parts. This diversity may not be apparent in current 
observations because the spread is contributed in part by a 
multitude of low-surface brightness features. Ongoing deep 
observations in the M3 1 should be able to test this prediction. 

- We have investigated the origin of high surface debris 
in different regions of the halo. The debris in the outer 
R > 10 kpc halo typically originates in single massive satel- 
lites. This result holds regardless of whether the debris is 
kinematically cold or if it is diffuse in phase-space. The outer 
regions give the best opportunity of reconstructing the for- 
mation history of stellar halos, as their stellar populations 
retain information about the nature of the mergers in a non- 
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degenerate form. 

- The origin of bright debris in the inner 10 kpc halo is 
more varied. Here tidal debris from a single massive pro- 
genitor or superposition of debris from several large or small 
satellites can result in similar surface brightness levels. Dis- 
entangling the complex formation histories of these fields is 
possible, provided that the interpretation of observational data 
takes into account the predictions of hierarchical models. 

- The brightest tidal debris (i.e. that can be observed with 
current surface brightness limits) still orbiting around M3 1 or 
Milky Way-type galaxies should be metal-rich. Therefore it 
should not be expected to resemble the characteristically old 
and metal-poor stellar populations of dwarf spheroidal satel- 
lites, nor those of the smooth underlying halo. 

- We illustrated the case of GSS in M31, whose progenitor 
was a massive dwarf galaxy accreted more recently than the 
main progenitors of the halo. The presence of intermediate- 
age stars and the absence of young stars suggest that the satel- 
lite could have been accreted as long ago as 6 - 7 Gyr If the 
satellite was accreted more than 1-2 Gyr ago, a continua- 
tion of the stream should exist in the halo of M31, about 1-2 
mag/arcsec^ fainter than the currently observed one. 
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